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Summary: Lattice modulations in the high magnetic field phase and 
close to impurities in spin-Peierls systems are considered and compared to 
experiment. Necessary extensions of existing theories are proposed. The 
influence of zero-point fluctuations on magnetic amplitudes is shown. 

The discovery of the first inorganic spin-Peierls substance CuGe03 by Hase et 
al. [H renewed the interest in the phases of spin-Peierls substances (for reviews, 
see pi). The constituting element of these substances are quasi one-dimensional 
spin chains coupled to the lattice degrees of freedom. The incommensurably 
modulated (I) phase is particularly interesting since this phase is the most com- 
plex one from the theoretical and from the experimental point of view. Already 
20 years ago much work was devoted to it (e.g. [3-9] ). Yet, detailed experimental 
investigations of the nature of this phase were not possible at that time. This 
has changed for CuGeC>3. 

X-ray experiments in the I phase permitted to detect the incommensurability 



of the lattice in fc-space 1 10 . The structure of the lattice modulation was inves- 
tigated by measuring the intensity of the third harmonic fl(i|| . The distribution 
of local magnetizations in CuGeC>3 was measured by NMR [fllj . Recently, the 
shape and the amplitude of the magnetic part of a soliton were deduced jn],[l|. 

Here we will review several aspects of the theoretical description of the I phase. 
Special emphasis is put on the discrepancies between the so far existing theories 
and the experimental findings and to resolving these discrepancies jl3| . 

The term "soliton" will be used for the combination of a zero in the modulated 
distortion and the concomitant localized spinon j|, The distortive soliton 
width £d is the width of the kink-like modulation. The magnetic soliton width £ m 
is the spatial width of the local magnetizations [l4| . The incommensurate 

modulation in the I phase is viewed as an equidistant array (lattice) of solitons. 

The schematic phase diagram in Fig. | shall serve as a guide for the following 
discussion. For the sake of concreteness we consider the hamiltonian 

H = J [(! + W+xS; + aSi+xSi-x] + y 8\ (0.1) 



which neglects the phonon dynamics (adiabatic treatment) since the distortions 
bi are numbers of which the average is zero. The nearest neighbour coupling J is 
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Figure 1 Generic phase diagram of a spin-Peierls system. U: uniform phase; 
D: dimerized phase; I: incommensurate phase. Solid line: second order transition; 
dashed line: weak first or second order. 

modulated by the <5,. Frustrating next-nearest neighbour coupling a J is also con- 
sidered. The spring constant K measures the elastic energy of a given distortion. 
Without magnetic field the spin chains are very susceptible towards dimeriza- 
tion Si = (— iyS. The energy lowering due to dimerization is anomalously large 
AE oc — <5 4 / 3 overcompensating the energy loss in elastic energy oc S 2 Jl6| . 

Applying a magnetic field to the dimerized (D) phase of a spin-Peierls system 
(T small) leaves the D phase unchanged until the critical field H c is reached. 
The stability stems from the singlet-triplet gap Atrip ■ Hence an upper bound for 
H c is g^H c = A^ip. In fact, the critical field H c is smaller. Its Zeeman energy 
represents only 80% of the gap A tr i P because the modulation also changes above 
the critical field. The system enters the incommensurably modulated phase. 



1 Sinusoidal Modulation 

Elastic X-ray scattering confirmed that the distortion in the I phase is in- 
commensurate. The deviation from dimerization \q — ir\ increases with increasing 
magnetic field. This can be illustrated in the unfrustrated XY model which corre- 
sponds via the Jordan- Wigner transformation to free fermions. The susceptibility 
towards distortion becomes maximum at q — 2kp since q = 2kp allows to create 
particle-hole pairs of vanishing energy. A distortion with q = 2fcp = 27rm + 7r is 
formed where m is the magnetization. Since this distortion couples the degener- 
ate states at — kp and at kp a gap appears there. Let us now gradually increase 
the interaction corresponding to the SfSf +1 terms at fixed particle number. For 
continuity no state can cross the gap such that the picture of a distortion at 2kp 
remains valid JIt], So we have 



\q — tt| = 27rm 



(1.2) 
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This relation is confirmed numerically p5| . 

Since experimentally it turned out that higher harmonics of the distortion are 



considerably suppressed it is plausible to start with (q given by (1.2) |l 



H = J^y. - *cos(gr i )]S i S <+1 . (1.3) 

i 

The distribution of the local magnetizations rm — (Sf) is found from NMR 
experiments [[Tl| |2). With some success, the experimental data were compared 
to a continuum theory |g. This theory, however, is based on a Hartree-Fock 
treatment where all Hartree and Fock terms are spatially constant. In Ref. |1£| 
it is shown that this is a too crude approximation reducing the physics to the one 
of a XY chain. The antiferromagnetic correlations found in this way are much 
smaller than those of an isotropic XYZ chain. This conclusion is corroborated 
by several works @, |20| But the spin isotropy of cuprates can hardly be 



questioned. To account for the smaller amplitudes it is proposed jl8|, |13| that 
experimentally only an effective magnetization mf is seen which is an average 

mf = (1 - 2j)ith + j(mi-i + m i+ i) . (1.4) 

The results for 7 = 0.2 agree well with experiment |lq] . The microscopic origin 
of the average is discussed in Sect. 4. 

The reason for strong local magnetizations around the zeros of the modulation 
(cf. Fig. p|) is found in the localization of a spinon. Each zero binds exactly one 
spinon 113]. Summing the m, around a magnetization maximum yields 1/2. 

The order of the transition D — » I can be determined by investigating the 
ground state energy E(m) as function of the average magnetization m |15| ], 
By means of a Legendre transformation E(h) — E(m) — hm, one obtains the 
dependence of ground state energy E{h) on the magnetic field h = gfi^H. It 
is found that a discontinuous jump for m — > occurs. This implies that the 
transition D — > I for fixed sinusoidal modulation is of first order. The mean 
square of cos(qri) jumps discontinuously from 1 to 1/2 if q deviates infinitesimaly 
from 7r [ p"5| since the ri are summed over integer values only. 

Experimentally, however, the observed first order jumps are much lower than 
those found for fixed sinusoidal modulation |^l|, |Io|, [n], g2|, |23| . 



2 Adaptive Modulation 

Since it was stated above that sinusoidal modulation alone does not account for 
the weak first order D — > I transition we turn to the full minimization of the 
ground state energy of d0.l|). Derivation with respect to Si yields 



= (S i+1 S i )-{{S j+1 S j ))+KS i , 



(2.5) 
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Figure 2 Symbols: DMRG-resuli_at K = 18J, a = 0.35. Upper panel: local 
distortions; solid line: from Eq. with 5 = 014, k d = 0.959, £ d = 10.5. 

Lower panel: local magnetizations; solid line: Eq. (2.6) with W = 0.21, R — 5.0, 
km = 0.992, £ m = 7.9. 



where ((•}} stands for the expectation value and the average along the chain. 
The double-bracketed term accounts for the constraint of the vanishing average 
of the Si. The minimization is done iteratively 

The generic result is depicted in Fig. |2[ The local magnetizations do not dis- 
play major differences to the results for sinusoidal modulation in ||l8) because 
the spinon localization is in essence determined only by the slope with which the 
modulation vanishes. Note that the envelope of m, is proportional to the proba- 
bility of finding a spinon at that site . Hence, the magnetic part of the soliton 
displays localization as for sinusoidal modulation. For the distortions the rela- 
tion (2.5) implies that the deviations from constantly alternating dimerization 
are also localized. 

The distortion belonging to an isolated soliton is a kink, i.e. the distortion 
between two solitons resembles the one in the D phase. This implies a crucial 
advantage over the sinusoidal modulation. For kink-like solitons the reduction of 
the mean square distortion is proportional to the soliton number. Hence, a low 
soliton concentration leads only to a small change of the energy such that E(m) 
is continuous (in the sense of Lipschitz) oam->0 Jig ]. 

The investigations in Ref. |15j of the model (0.1) yielded a continuous phase 
transition even though the magnetization grows very quickly above the critical 
field m oc — l/hx{H — H c ). Most of the results of the continuum theories comply 
also with a phase transition of second order || ||, [?], |) . Solely Horovitz mentions 
the possibility of soliton attraction in an early work || implying a first order 
transition. Buzdin et al. expect a first order transition at T > [Q. 

In fact, the details of the models matter. Cross J23 argued already that an 
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elastic energy with dispersion K (q) being minimum at q = tt leads to a first 
order transition. The positive curvature of K(q) around q = n suppresses higher 
harmonics in the distortion. Hence, sinusoidal modulation is favoured. The con- 
comitant concavity in E(m) at low magnetization m implies phase separation via 
the Maxwell construction. So the phase transition is first order This finding 
is in accordance with the conclusion from a phcnomenological Ginzburg-Landau 
description |p5| that the D — > I transition is generically of first order. 

At the transition the distance between two solitons is rather large so that the 
difference between sinusoidal and adaptive modulation matters most. At higher 
soliton concentrations the adaptive modulation becomes more and more sinu- 
soidal but with a concentration dependent amplitude. The mean square distor- 
tion could be determined from the elastic lattice constants. These experimental 
results agree well with the predictions based on the model (p.l[), see p6 . 

The continuum theories applying to the isotropic Heisenberg chain 27) pro- 
vide the following results (details in Ref. (l3|; sn, cn, dn: elliptic Jacobi functions) 
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The fits in Fig. | are based on Eqs. (|2.6| 



n/dKm/d 



if). Identity (2 



-,fcm)J (2.6) 

(2.7) 

(2.8) 
(2.9) 

(2.10) 

is not complied 



with, see Sect. 3. Otherwise no agreement would be obtained. Relation (2.9) is 
imposed on the fits whereas Eq. ( |2.1C ) serves as check. It is fulfilled within 4%. 
The fact that £d/£m ~ 1-33 is considerably above unity complies nicely with the 
experimental findings. Elastic X-ray scattering [jl0| found £d = 13.6 ± 0.3 while 
the NMR investigations provide £ m « 10 |H| close to the transition. 

There is also another way to introduce solitons in a spin-Peierls system than 
the application of magnetic field. Doping non-magnetic impurities in a spin- 
Peierls systems cuts the infinite chains into finite chain segments p8| . In a num- 
ber of works p9| it has been shown that each impurity frees one spinon which 
is situated either before or after the impurity on the chain. Assuming a fixed 
dimerization it is easy to see that the spinon is bound to its generating impurity 
Q in accordance with experimental results ||l| . 

But in a spin-Peierls system the change of the modulation has to be taken 
into account, too. This is done by introducing 



H 



i>0 



K 

[1 + 5i)SiSi+i + aSjS i+ 2 + — 5 t + f5i(—iyShuik 



(2.11) 
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Figure 3 Local distortions (enhanced) and magnetizations in the vicinity of an 
impurity (chain end) at site -1 for K = 18.064, a = 0.35 and / as displayed. 



such that the impurity is at site -1. The important amendment compared to H 
in Eq. ( p.l| ) is the last term. If the spinon moves away from the impurity the dis- 
tortion pattern is changed between impurity and spinon. Due to an elastic inter- 
chain interaction (parametrized by /) a coherent distortion pattern throughout 
the whole three-dimensional system is preferred. A deviation from this p attern 
is energetically unfavourable. So one is led to include the last term in Eq. ( 2.11 ) 
assuming that the adjacent chains are dimerized as in the unperturbed D phase. 
The relation K = Kq + f ensures the consistency of the distortion amplitude 
with its bulk value. 

Fig. ^ displays the results for various elastic interchain interactions. The soliton 
is not localized at the impurity but at a certain distance. For lower / values the 
soliton resembles very much the ones in the I phase (cf. Fig. ||). On increasing / 
the soliton is squeezed more and more towards the impurity. Analogous results 
can be found by QMC §§, too. 

Sure enough, the soliton is bound to the impurity confirming previous ideas 
|p8| , [n]. Moreover, the first excitations for the same distortions and in the same 
spin sector are found below the singlet-triplet gap at 52% of A tl i P for / = 0.01 
and at 64% for / = 1. This matters for the spectroscopic analysis. 



3 Local Renormalization 

Most remarkable is the difference between the distortive soliton width £d and the 
magnetic soliton width £ m in the numerical results (cf. Fig. ||). It amounts to 30% 
at a = 0.35 depending mainly on the frustration for low soliton concentrations 
. The challenge is to extend the existing continuum description to account for 
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Figure 4 Local fluctuations as function of x. The deviation Act from the ground 
state value (solid line) and the renormahzing factor (dashed) is shown. 

this fact. Let us revisit the semiclassical treatment of the bosonized description 
of the spin-Peierls problem Q . Minimizing the total energy the variation of the 
distortion 8{x) leads to 

S(x) Ke- 2CT cos(20 class ) (3.12) 

where u := ((f) 2 ) denotes the renormahzing fluctuations of the local bosonic field 
4> about the classical field c iass- A soliton corresponds to a solution where <^ c iass 
increases by ir in a kink- like fashion. If a is assumed to take the spatially constant 
value that it has in the ground state one obtains 

S(x)/S = cos(20 c i ass ) = tarm(x/£) , (3.13) 

wherein £ is given by the ratio ^s/A t ri P of the spin wave velocity and the gap. The 
alternating component of the magnetization a(x) is proportional to sin(2</> c i ass ). 

Hence, one has a(x) oc yl — tanh 2 (x/£) = 1/ cosh(a;/£). 

But the presence of the soliton induces a deviation Aer from its ground state 
value. Fig. |J displays a generic result for this deviation. It is calculated on top 
of the solution of Nakano and Fukuyama pf . The alternating component 



a(x) oc y 1 - exp(4Aa) tanh 2 (:r/£) (3.14) 

is thus indeed narrower than before due to the spatial dependence of the renor- 
malization factor. The result in Fig. ^ is only a first step since the influence of 
the altered magnetic behaviour is not included. Yet it is clear that the renormal- 
ization is a local quantity and that this is the origin of the difference between £d 
and Cm- 



8 



Uhrig 



4 Phasons 



In Sect. 1 the discrepancy between theoretical and experimental amplitudes of 
the alternating local magnetizations was solved by an averaging procedure ( |l.4|) . 
Passing to adaptive modulations does not change the amplitudes much ]l3fl . 
Hence, one still has to find the microscopic origin for the averaging. 

Non-adiabatic effects are in fact responsible. The soliton lattice oscillates 
about its equilibrium positions. These oscillations are best understood in a con- 
tinuum description which is well justified if the typical length £ is noticeably 
larger than the lattice constant. In such a continuum description the modula- 
tion 5(x) can be shifted along the chains without energy cost. This continuous 
translational invariance is spontaneously broken by the soliton lattice. Hence, 
there are massless Goldstone modes, the so-called phasons. They are analogous 
to the phonons of a crystal lattice except that they do not have three branches 
but only with one. While the atoms of a crystal lattice can be shifted in all three 
spatial directions the solitons can be shifted only along the chains. 

Albeit the ideal spin-Peierls system is magnetically one-dimensional the dis- 
tortions on different chains are elastically coupled. Thus, the phasons are gov- 
erned by a 3D, though anisotropic, dispersion. The dispersion parameters are 
determined from the anisotropy of the correlation lengths assuming a Ginzburg- 
Landau description f2q| . The corresponding T 3 term in the specific heat has been 
measured and theory and experiment agree astonishingly well f25fl . 

The zero point motion and the excited motion of phasons lead to the aver- 
aging (-L4). Let us denote the adiabatic result for the local magnetizations by 
rrii — a(ri) cos(irri) + u(ri) where is the component along the chains, a{r{) the 
alternating component and u(ri) the uniform component of the magnetizations. 
A local shift can be implemented by replacing -kti — > -kti + 0(r,) where 0(i"j) 
denotes a phase shift operator. On the long time scales of a NMR measurement 
one measures 

m^ xp = {mi) = 0(^)7' cos(7rr 4 ) + u(r;) (4-15) 

with 7' := exp ^— ^ ^ i (9 2 (r i ))^ < 1. So there is an amplitude reduction 

engendered by the local fluctuations. The factor 7' is comparable to a Debye- 
Waller factor. 

Using the values fixed previously §| yields 7' = 0.16 cxp(-(T/T*) 2 /2) with 
T* k 16. 9K Jl3| . Eq. (1.4) is retrieved by estimating a(r) and u(r) from the 
discrete values m, by alji) — m,i/2 — (rrii^i + mj + i) /A a nd by u{ri) = m j/2 + 
(rrii-i + to;+i)/4. Inserting these formulae into Eq. ( |4.15| ) yields Eq. ( |l.4| ) with 
7 = (1 — 7')/4- At T = 0, 7 takes the value 0.21 in accordance with experiment 

ram. 
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5 Conclusions 

In this report the modulated phases of spin-Peierls systems were discussed. Such 
modulations are induced either by magnetic field or by impurities. In both ways 
the singlet pairing in the D phase is broken and spinons are freed. The lattice 
distortion adapts to the spinon by forming a zero to which the spinon is bound. 
This new entity constitutes the spin-Peierls soliton. 

The order of the transition D — > I phase on increasing field depends on model 
details. Imposed sinusoidal modulation leads to a pronounced first order transi- 
tion. Allowing the system to choose an optimum modulation makes the transition 
continuous if the elastic energy is wave vector independent. If the elastic energy 
itself pins the modulation to tt a weak first order transition is found. 

Doping induced solitons are bound to their generating impurity. The distortion 
pattern between impurity and soliton is not coherent with the bulk pattern. 
This costs energy which acts as a confining potential. Binding occurs for which 
experimental evidence exists |3lJ . 

The difference between magnetic and distortive soliton width could be traced 
back to the so far neglected spatial dependence of the renormalizing local fluc- 
tuations. A fully self consistent analysis is in progress. 

The reduction of the alternating magnetic amplitude due to phasons provides 
striking evidence for the importance of the lattice dynamics. The inclusion of 
non-adiabatic effects on top of an otherwise adiabatic calculation might still 
be unsatisfactory. So other approaches to non-adiabatic behaviour should be 
extended to the I phase |34| . 

I thank C. Berthier, J. P. Boucher, B. Biichner, T. Lorcnz, M. Horvatic, Th. 
Nattermann for fruitful discussions, F. Schonfcld for reliable numerical work, E. 
Miillcr-Hartmann for generous support and G. Giintherodt's group for intense 
collaboration. Financial support of the DFG by the SFB 341 is acknowledged. 
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